Superconductors of mixed order parameter symmetry in a Zeeman magnetic field 
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We study the effect of Zeeman magnetic field in the superconducting phase with two component 
order parameter scenario, such as, d x 2_ y 2 + e 1 a, where a — d xy , s. This scenario is equivalent to 
applying magnetic field parrllel to Cu — O planes. In a weak magnetic field, which does not cause 
much change to the predominant d-wave, supresses the minor a component leading to pure d-wave 
phase. This observation is in contrary to the effect of magnetic field applied in the c-direction to 
the Cu — O planes which is believed to induce a minor component a to d-wave superconductors. We 
also show that the response of such superconductors to a weak Zeeman magnetic field can be quite 
different depending on the phase 8 of the minor component (a). 
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I. INTRODUCTION 



Although the nature of the superconducting pair wave 
function in high -T c cuprates is not yet known strong ev- 
idences of a major d x 2_ y 2 symmetry exists |j]-|3|]. Experi- 
ments sensitive to the internal phase structure of the pair 
wave function reported a sign reversal of the order pa- 
rameter supporting d wave symmetry M. Most recently 
from various experiments and theory it appears that the 
pairing symmetry of these family could be a mixed one 
like d x 2_ y 2 + e l9 a where a could be something in the s 
wave family or d xy . There were early questions from tun- 
neling experiments regarding the pure d-wave symmetry 
H as the data supports an admixture of d and s-wave 
components due to orthorhombicity in YBCO Pos- 
sibility of a minor but finite id xy symmetry alongwith 
the predominant d x i_ y i has also been suggested []1 0f| in 
connection with magnetic defects or small fractions of 
a flux quantum <i>o = hc/2e in YBCO powders. Similar 
proposals came from various other authors in the context 
of magnetic field, magnetic impurity, interface effect etc. 
|jll"|-|l3f The experimental result by Krishana et al., was 
interpreted as a signature of induction of a minor com- 
ponent eg, id xy or is in a d-wave superconducture with 
the application of magnetic field along the c axis. 

In this work, we study in details the effect of a weak 
Zeeman magnetic field on superconductors with mixed 
order parameter symmetry like A(fe) = 2 _ 2 + e i6 a 
with a = d xy , s for arbitrary 8. It is well known that such 
superconductors with 8^0 and a ^ corresponds to 
broken time reversal states (BTRS). These BTRS states 
lift the directional degeneracy of charge currents by ad- 
mixing a subdominant cv-wave component to the d-wave 
pairing state and a spontaneous finite current appears 
jl4j]. An application of Zeeman magnetic field can lift 
the spin degeneracy leading to suppression of BTRS ; a 
pure d-wave occurs with increasing magnetic field. For 
mixed symmetry as above with 8 — that preserves time 
reversal symmetry and are nodeful respond differently to 



the Zeeman field as compared to 8 ^ states. For node- 
full 8 = state, the local gap A(fc) of small magnitude 
over the Fermi surface may be destructed with the appli- 
cation of Zeeman field leading to a paramagnet pocket. 
This although true for 8 ^ states, but such states cor- 
respond to fully gapped situation all over the Fermi sur- 
face which causes weak response to the magnetic field. A 
clear picture on the above will be demonstrated in this 
article. It may be mentioned that the high temperature 
superconductors are quasi two dimensional in nature and 
therefore, a magnetic field parallel to the Cu — O plane 
does not couple to the orbital motion of the electrons 
in the plane. Therefore, we shall not consider spin-orbit 
interaction in this work. 

In connection with the discussion of order parameter 
symmetry in cuprates, we would further like to mention 
that the proposal of mixed order parameter symmetry 
got the correct momentum when experimental data on 
longitudinal thermal conductivity by Krishana et al, 
of Bi2Sr2CaCu208 compounds and that by Movshovich 
et al, showed supportive indication to such propos- 
als. There are experimental results related to interface 
effects as well as in the bulk that indicates mixed pairing 
symmetry (with dominant d-wave) flTjjfl , thus providing 
a strong threat to the pure d wave models. There were 
early orginal works as regards to the modifaction of su- 
perconductivity due to application of Zeeman magnetic 
|L5[ and very recently, the Zeeman suppression was dis- 



cussed in mesoscopic systems |16 



II. MODEL CALCULATION 

The free energy of a two dimensional planar supercon- 
ductor with arbritary pairing symmetry in presence of a 
magnetic field may be written as, 
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where E% = y/(e k - A*) 2 + | A fc | 2 + a(gfi B /2)B are the 
energy eigen values of a Hamiltonian that describes su- 
perconductivity, is the magnetic moment of the 
electrons. This includes the assumption that the Zeeman 
field raises/lowers the energy of the spin up/down quasi- 
particle states. We minimize the free energy, Eq. ^ i.e, 
dF/d | A | = 0, to get the gap equation as, 
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where e k is the dispersion relation taken from the ARPES 
data j^] and p the chemical potential will control band 
filling through a number conserving equation given be- 
low. 

Since the applied Zeeman field modifies the SC quasi- 
particles of spin up and down differently, their occupation 
probablities are also modified. The number conserving 
equation that controls the band filling through chemical 
potential, /i in presence of Zeeman field is given by, 
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where h = (g/j,B/2)B. Let us consider that the overlap 
of orbitals in different unit cells is small compared to the 
diagonal overlap. Then in the spirit of tight binding lat- 
tice description, the matrix element of the pair potential 
used in the SC gap equation Eq. (g) may be obtained as, 

V(q) = E Vf 1 ^ = V ° + Vif d (k)f d (k') + V ig {k)g{k') 
5 

+V 2 f d ^(k)f d ^(k') + V 2 f s * v (k)f s * v (k') (4) 



where in the first result of the equation (||) Rs lo- 
cates nearest neighbour and further neighbours, <5 la- 
bels and V n , n = 1,2 represents strength of attraction 
between the respective neighbour interaction. The first 
term in the above equation Vq refers to the on-site in- 
teraction which has an effective attractive value giving 
rise isotropic s wave. The second and third terms are 
responsible for d and extended s wave symmetry super- 
conductivity whereas the A th and the 5 th terms are re- 
sponsible for d xy and s xy symmetries respectively. We 
restrict only to singlet pairing states (i.e, A(k) = A(— k)) 
as applicable for high temperature superconductors. The 
momentum form factors are obtained as, 

f d (k) = cos(k x a) — cos(k y a) 

g(k) = cos(k x a) + cos(k y a) 

f d *y(k) = 2 sin(fc x a) sm(k y a) 

f Sxv (k) = 2cos(k x a) cos(k y a) (5) 

For two component order parameter symmetries as men- 
tioned above, we substitute the required form of the po- 
tential and the corresponding gap structure into the ei- 
ther side of Eq. (0) which gives us an identity equation. 



Then separating the real and imaginary parts together 
with comparing the momentum dependences on either 
side of it we get gap equations for the amplitudes in dif- 
ferent channels as, 
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where j = 1,2 corresponding to two components d and 
s or d and d xy symmetries. Considering mixed symme- 
try of the form A(k) = Ad x2 _ y2 (0)/ d (fc) + e w A s {0) one 
identifies A, = A^^O)* A 2 = A s (0), ft = f d (k), 
fl = 1 and V 1 = Vl, V 2 = V in Eq. (4). Similarly, for 
mixed symmetries of the form A(k) — Ad 2 _ 2 (0)f d (k) + 

e* e A d ^{Q)f d *y A 2 = A dxH (0), /<? = and V 1 = V u 
V 2 = V? of Eq.(4). The potential required to get such 
pairing symmetries are discussed in Eq. (Q). 

We solve self-consistently the above three equations 
(Eq.^| and Eq.||) in order to study the phase diagram of 
a mixed order parameter superconducting phase in pres- 
ence of Zeeman magnetic field. The numerical results 
obtained for the gap amplitudes through Eqs. (^S) will 
be compared with free energy minimizations via Eq. (|l|) 
to get the phase diagrams. 



III. RESULTS AND DISCUSSIONS 



We present in this section our numerical results for a 
set of fixed parameters, e.g, a cut-off energy i7 c = 500 
K around the Fermi level above which superconducting 
condensate does not exist, a fixed transition tempera- 
ture of the minor component T"(h — 0) = 24 K and 
the bulk T c — 85 K determined by the d-wave order 
parameter. In figures 1 and 2 we present results for 
AO) = A d ^ 2 y2 (0)/ d (/c) + e ie a symmetries for 6 = tt/2 
and 9 = respectively. Such symmetries would arise 
from a combination of two component pair potentials 
(2 nd , 1 st ), (2 nd , 4 th ) terms in Eq. | for a = s k d xy 
respectively. The amplitudes of extended s-wave states 
like S2.2_i.y2, s xy are found to be finite only towards very 
low band filling, p ~ for 9 = tt/2 hence does not cause 
any mixing with the predominant ef-wave. Therefore, we 
shall discuss only the results of 9 = and 9 = tt/2 for 



a = s & d xy . These two phases of 9 can cause important 
differences. It is known that for any 9^0, time reversal 
symmetry is locally broken Jig] at lower temperatures 
with the onset of the secondary component which cor- 
respond to a phase transition to an fully gapped phase 
for h = 0, from a partially ungapped phase of d x i_ y i 
symmetry. On the other hand, the 9 = phase still re- 
mains nodeful, although the nodal lines shifts a lot from 
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the usual k x — k y lines of the d x 2_ y 2. In Fig. 1(a) and 
(b) we present the temperature dependencies of the order 
parameters in the complex mixed symmetry. Curves cor- 
responding to the (i-wave channels are represented with 
thinner joining lines of different styles whereas the minor 
s component is denoted through that of same style but 
with thicker lines (same strategy will be carried out in 
other figures as well). For a zero Zeeman field (h = 0), 
the amplitude of the d-wave is suppressed with the onset 
of the minor s component at T = 24 K which leads to a 
kink like structure (cf. the thin solid curve in Figl (a)). 
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FIG. 1. (a) Amplitudes of the Ad 2 _ 2 and A s at a fixed 
band filling p = 0.85 as a function of temperature (in K) for 
9 — 7r/2 (i.e, d x 2_ y 2 + is) phase in various values of Zeeman 
field (h). With the application of weak magnetic field while 
the d x 2_ y 2 remains almost unaffected the transition temper- 
ature (T C Q ) of the minor a = s component is suppressed 
strongly. At a field value of h — 0.014 eV the s-wave shows 
a first order transition, (b) For field values h > 0.016 eV the 
minor s-channel is completely suppressed leading to a pure d 
wave order parameter. The thermal behaviour of d-wave su- 
perconductivity in presence of the Zeeman field is presented 
in Fig. 1(b). While the amplitudes at T = K remains almost 
unaffected with field, the transition temperatures get affected. 
A magnetic field induced first order transition is observed at 
h = 0.04 eV. 

With application of weak Zeeman field the transition 
temperature of the minor s state decreases while the zero 



temperature magnitude remains the same. This causes 
a shift in the kink like structure in the thermal depen- 
dence of the d wave channel towards lower temperature 
with increasing field and hence a small enhancement in 
the d-wave with field at lower temperature occurs. This 
point will be clearer from Fig. 5(b) as discussed latter. 
At a field value h c = 0.016 eV, the s wave component 
is completely suppressed leading to a pure d wave phase. 
Thus we have a magnetic field induced transition at lower 
temperature from fully gapped phase of the d + is state 
to a partially gapped phase of the <i-wave. Therefore, 
in absence or very low magnetic field, there is a transi- 
tion from a partially gapped phase of the <i-wave at higher 
temperature to a fully gapped d + is phase at lower tem- 
peratures. With increasing field the fully gapped phase 
region with respect to temperature decreases and brings 
back the ungapped phase of the d-wave. These phase 
transitions will have important bearings in the thermo- 
dynamic and transport properties. 
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FIG. 2. Same as that of figure 1(a) except 9 — (i.e, 
d x 2_ y 2 + s symmetry). For the d-wave channel, its ampli- 
tude is found to increase below T c s in contrast to 9 — n/2 
case where the d-wave amplitude is suppressed. The s-wwave 
channel suffers drastic suprresion in T c s as well as the zero tem- 
perature amplitude in contrast to that in Fig. 1(a). Note, the 
critical field (h c ) at which the s component is completely sup- 
pressed is 0.013 eV in contrast to 0.016 eV in case of 9 = 7r/2 
(cf. Figla). 

In Fig. 1(b) we present the paramagnetic state of the 
d-wave. This phase has also been qualitatively investi- 
gated recently by Yang and Sondhi with possibility 
of pairing with finite momentum. We therefore restrict 
to present only the details thermal dependence of the d- 
wave superconductivity in presence of Zeeman magnetic 
field that were not discussed. We show that with increas- 
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ing field (at lower fields) the zero temperature magnitude 
of the d-wave remains unchanged whereas the T c is re- 
duced. At higher field e.g, h — 0.04 one sees a first order 
transition from superconducting state to the normal state 
with respect to temperature. This behaviour causes a 
magnetic field induced enhancement of the 2A//csT c ra- 
tio. This ratio is crucial for many physical properties 
like specific heat jump etc. and hence expected to have 
drastic effect with magnetic field. 
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FIG. 3. Amplitudes of the Ad 2 _ 2 an d Ad xy at a fixed 
band filling p — 0.85 as a function of temperature (in K) for 
9 — it/2 (i.e, d x 2_ y 2 + id xy in various values of Zeeman field 
(h). With the application of weak magnetic field while the 
d x 2_ y 2 remains practically unaffected the transition temper- 
ature (T™) of the minor a = d xy component is suppressed 
strongly. Similar to that of Fig. la, a first order transition is 
seen in the minor d xy channel for h = 0.014 eV. 

In Fig. 2 we describe the thermal behaviour of super- 
conductors with d + s symmetry in presence of magnetic 
field. First of all, in absence of magnetic field, the ther- 
mal behaviors at lower temperatures is quite different (as 
we have seen in Fig. 1), the s-wave gap opens very fast 
below T = 24 K and also induces a growth to the d at a 
faster rate than that above T — 24 K. While the s-wave 
has about three times the zero temperature value com- 
pared to d + is phase, the d wave also have quite larger 
value. With application of a very small magnetic field 
the s-component is suppressed largely; both its T c and 
the zero temperature gap being suppressed. The d-wave 
gap magnitude is also suppressed while its T c remained 
practically unchanged with such small values of the mag- 
netic field. More importantly, although both the d and 
the s channel has larger magnitude in the d+s phase the 
critical field (h c — 0.013 eV) at which the s component 
vanishes completely is smaller compared to that for the 



6 = 7r/2 phase (h c = 0.016 cV) of the mixed symme- 
try Distinctly, the response of the Zeeman field to the 
d + s superconductors is more pronounced than the d + is 
superconductors. This study therefore also revealed the 
importance of the phase of the minor component in d- 
wave superconductors with a pratical example of effect 
of Zeeman magnetic field, for the first time. In order to 
establish the importance of the phase of the minor com- 
ponent we also study the effect of Zeeman magnetic field 
in d x 2_ y 2 +id xy and d x 2_ y 2 + d xy symmetry superconduc- 
tors respectively. The qualitative as well as quantitative 
behaviors remain almost same as that shown in figures 1 
and 2. Thus establishing that the response of the super- 
conducors with mixed order parameter symmetries with 
9 = phase of the minor component is stronger. 
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FIG. 4. Same as that in figure 3 except for 9 = i.e 
d x 2_ y 2 + d xy phase that preserves the time reversal symme- 
try. The notable difference is that the minor d xy component 
although have a magnitude at T = K and h — three times 
larger than that for — -k/2 case, the minor component is 
suppressed at a lower critical field h c = 0.013 eV. 



In figures 5 and 6 we concentrate on studying behavior 
of such mixed order parameter symmetry superconduc- 
tors in presence of magnetic field as a function of band 
filling (p) for d xy and s as minor components respectively. 
With application of a small field the minor component 
a = d xy or s is suppressed only at lower fillings and then 
with increasing field it is suppressed in the optimal dop- 
ing regime suddenly Therefore, with increasing field one 
finds mixed symmetry region as well as pure d-wave re- 
gion with respect to filling - a non-uniform superconduc- 
tivity. This behavior depends on the nature of the minor 
component. For a = s, the mixed symmetry is possible 
around half-filling and at around p = 0.7 in an interme- 
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diate field value. With increasing field mixed symmetry 
regions around both the band fillings shrinks and at a 
field value of h = 0.03 the s-wave vanishes leading to a 
pure (i-wave phase. It may be noticed the (i-wave channel 
is enhanced (cf. Fig 5(b)) in the optimal doping region 
with intermediate field values as was also mentioned ear- 
lier while discussing Figl (a). This however does not oc- 
cur in case of d xy . For d xy minor component, the minor 
component is suppressed strongly only from the lower 
filling with increasing field and at an intermidiate field 
the mixing is possible only near half filling. The (i-wave 
boundary towards lower fillings as well as around half- 
filling also shrinks with increasing field values in contrast 
to a = s. One always see a sharp transition from a mixed 
phase to a pure (i-wave phase or otherway around with 
respect to band filling irrespective of the minor compo- 
nent in a given magnetic field. At larger fields (h > 0.03) 
when the minor component is suppressed completely i. e, 
one has a pure (i-wave phase, the (i-wave phase also show 
sharp transition from (i-wave superconducting state to a 
normal state (cf. Fig. 6). 
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In Fig. 6 we study the same as that in Fig. 5 for d+ s 
and d+d xy symmetries. In contrast to the 9 = tt/2 phase 
of the minor component (as in Fig. 5), the minor compo- 
nents have very large values in the 9 = phase although 
have the same T" as mentioned earlier, in absence of the 
magnetic field. With very small magnetic field such large 
amplitudes of the condensation in the minor channel is 
strongly suppressed (see Fig 6(b) specially). The nature 
of suppression, as in Fig. 5, is different for different minor 
component. For example, d xy is suppressed only from the 
lower filling whereas the s-wave is suppressed both from 
the lower filling as well as around the optimal doping. 
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FIG. 6. Same as figure 5 except 8 = 0. Stronger and rapid 
suppression of the minor channel with field in contrast to that 
in Fig. 3. A visible change in the dominant d-wave channel is 
also seen in contrast to that in Fig. 3. Note, in case of d + s 
picture, (b), the d-wave channel shows plateau with respect 
to band filling. 



FIG. 5. Zero temperature amplitudes of the different pair- 
ing channels as a function of band fillings in the d x 2_ y 2 +e 10 a 
picture for several fields with 9 = n/2. An increase in the 
amplitude of the d-wave channel with some particular field 
values in the d x 2_ y 2 +is picture near optimal doping is worth 
noticing. With field increasing the order parameter symmetry 
changes to pure d-wave at optimal and larger doping. 



Overall, it is very distinct that the magentic field af- 
fects the 9 = mixed phase more stongly than the 
9 = tt/2 phase. This is due to the fact that the response 
of applied Zeeman field is paramagnetic with destruc- 
tion of superconductivity over parts of the Fermi surface 
where the Zeeman field exceeds the local magnitude of 
the fc-dependent gap resulting a spin polarization in the 
normal electrons. For 9 = n/2 phase, the nodes are miss- 
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ing and the Fermi surface is gapped all over, although the 
gap will have local minima. Thus response of the 9 = n/2 
phase is weaker compared to the 8 = phase. Noticiably, 
the critical value at which the minor component vanishes 
completely in all band fillings for the 9 = phase is 
h c = 0.025 eV whereas that for 8 = ir/2 is h c = 0.03 eV. 
(The value of h c depends on p, as well as a and 8). 

Above h c = 0.025, the paramegnetic state of the d- 
wave superconductors are also presented. The region 
of (i-wave superconductivity shrinks with increasing field 
around the band filling p = 0.82. The change from d- 
wave to spin polarized normal state is very sharp with 
respect to the band filling. In the d + s state a plateau 
has been observed in the d channel in weak or zero field, 
similar to the behaviour known for the YBCO systems 
|fl9|| . It may be mentioned that Krishana et ah, found the 
longitudinal thermal conductivity of Bi2Sr2CaCu20s+5 
at lower temperatures (5K to 20 K) decreases with the 
increase in magnetic field applied along c-axis. Above 
a critical value of the magnetic field Hk(T), the ther- 
mal conductivity cease to change with the magnetic field 
and develops a plateau. It was proposed [jl2) that the 
d x 2_ y 2 pairing state is unstable against the formation of 
d+e i6 a (where a = s, d xy ) in presence of Hk(T) such that 
the loss of quasiparticle transport in the thermal conduc- 
tivity can be explained. In contrast, we started with a 
d + e l6 a picture and application of Zeeman field (which 
may be mapped as application of magnetic field parallel 
to the 2D Cu — O plane) without considering the orbital 
effect did not find any enhancement in the condensation 
of the minor channel but suppression leading to a pure 
paramegnetic d state. 

IV. SUMMARY 

In summary we have performed a detailed study on the 
effect of Zeeman magnetic field on mixed pairing sym- 
metry with predominant d-wave which seems to be very 
promising symmetry for the high T c systems. Thus we 
have described the paramagnetic state in the mixed sym- 
metry superconductors and subsequently in the d-wave 
superconductors. In particular we established that the 
phase of the minor component mixed with predominant 
d-wave is of immense importance. The 9 = phase mi- 
nor component symmetry responds to Zeeman field more 
profoundly than the 9 = n/2 of the minor component. 
It will be very interesting to calculte the specific heat, 
Magnetization, density of states as a function of mag- 
netic field using this model. We argued that the orbital 
effect is secondary when a magnetic field is applied par- 
allel to the conducting plane. This may indicate that 
the experimental observation by Krishana et al., involve 
strong coupling of spins to orbitals due to application of 
magnetic field perpendicular to the plane at lower tem- 
peratures. However, the order parameter alone does not 



completely determine the thermodynamic property of a 
system. In Ref. |l7j|, some estimations and scaling rela- 
tions in the change of various physical properties due to 
Zeeman field are given for a pure d-wave superconductor. 
It turns out, a weak Zeeman field does little to the order 
parameter but may profoundly affect the thermodynamic 
property of a pure d x 2_ y 2 supercondcutor. Such effects 
will presumably also remain for mixed superconductors 
with 9 = as the gapnodes similar to d x 2_ y 2 remains, 
its effect in the 9 = it/ 2 phase has to be studied more 
carefully. 



V. ACKNOWLEDGMENTS 

A large part of this work was carried out at the Insti- 
tute) de Fisica, Universidade Federal Fluminense, Brazil 
and was financially supported by the Brazilian funding 
agency FAPERJ, project no. E-26/150.925/96-BOLSA. 



[1] D. L. Cox and M. B. Maple, Phys. Today, 32, Feb. 1995. 

[2] D. J. Van Harlingen, Rev. Mod. Phys, 67, 515 (1995). 

[3] D. J. Scalapino, Phys. Rept. 250, 329 (1995). 

[4] D. A. Wollman, D. J. Van Harlingen, W. C. Lee, D. M. 
Ginsberg and A. J. Leggett, Phys. Rev. Lett. 71, 2134 
(1993) ; D. A. Brawner and H. R. Ott, Phys. Rev. B 50, 
6530 (1994) ; C. C. Tsuei, J. R. Kirtley, C. C. Chi, Lock 
See Yu-Jahnes, A. Gupta, T. Shaw, J. Z. Sn and M. B. 
Ketchen, Phys. Rev. Lett. 73, 593,(1994) ; W. N. Hardy, 
D. A. Bonn, D. C. Morgan, R. Liang and K. Zhang, 
Phys. Rev. Lett. 70, 3999 (1990); L. A. de Vaulchier, 
J. P. Vieren, Y. Guldner, N. Bontemps, R. Combescot, 
Y. Lema itre and J. C. Image, Europhys. Lett. 33, 153 
(1996). 

[5] Dong Ho Wu, J. Mao, S. N. Mao, J. L. Pong, X. X. Xi, 
T. Venkatesan, R. L. Greene and S. M. Anlage, Phys. 
Rev. Lett. 70, 85 (1993) ; S. M. Anlage, Dong-Ho Wu, 
J. Mao, S. N. Mao, X. X. Xi, T. Venkatasan, J. L. Peng 
and R. L. Greene, Phys. Rev. B 50, 523 (1994) ; A. F. 
Annett, N. Goldenfield and A. J. Leggett, Jr. Low Temp. 
Phys. 105, 473 (1996). 

[6] A.G. Sun, D. A. Gajensk, M. B. Maple and R. C. Dynes, 
Phys. Rev. Lett. 72, 2267, (1994). 

[7] M. B. Walker, Phys. Rev. B 53, 5835, (1996). 

[8] C. ODonovan, D. Branch, J. P. Carbotte and J. S. Pre- 
ston, Phys. Rev. B 51, 6588, (1995). 

[9] M. R. Norman, M. Randeria, H. Ding and J. C. Cam- 
puzano, Phys. Rev. B 52, 615, (1995). 
[10] M. Sigrist, D. B. Bailey, and R. B. Laughlin, Phys. Rev. 

Lett. 74, 3249, (1995). 
[11] R. B. Laughlin, Phys. Rev. Lett. 80, 5188, (1998). 
[12] K. Krishana, N. P. Ong, Q. Li, G. D. Gu, and N. 
Koshizuka, Science, 277, 83, (1997); K. A. Kouznetsov 



6 



et ai, Phys. Rev. Lett. 79, 3050 (1997) ; A. V. Bal- 
atsky, Phys. Rev. Lett. 80, 1972 (1998); R. Movshvich, 
M. A. Hubbard, M. B. Salamon, A. V. Balatsky, R. 
Yoshizaki, J. L. Sarrao and M. Jaime, Phys. Rev. Lett. 
80, 1968 (1998); R. B. Laughlin, P hys. Rev. Lett. 80 
5158 (1998) ; T. V. Ramakrishnan, |sond-mat/9803069| ; 
Haranath Ghosh, Europhys. Lett. 43, 707 (1998). 
[13] M. Covington, G. Aprili, E. Paraoanu, L. H. Greene, F. 
Xu, J. Zhu and C. A. Mirkin, Phys. Rev. Lett. 79,277 

(1997) ; M. Forgelstrom, D. Rainer and J. A. Sauls, Phys. 
Rev. Lett. 79, 281 (1997); M. E. Zhitomirsky and M. 
B. Walker Phys. Rev. Lett. 79, 1734 (1997); R. J. Kel- 
ley, C. Quitmann, M. Onellion, H. Berger, P. Almeras, 
G. Margaritondo, Science 271, 1255 (1996) ; J. Ma, C. 
Quitmann, R. J. Kelley, H. Berger, G. Margaritondo, M. 
Onellion, Science 267, 862 (1995). 

[14] For a review see M. Sigrist, Prog. Theor. Phys. 99, 899 

(1998) . 



[15] A. M. Clogston, Phys. Rev. Lett. 3, 266 (1962); B. S. 
Chandrasekhar, Appl. Phys. Lett. 1, 7 (1962); P. Fulde 
and R. A. Ferrell, Phys. Rev. 135, A550 (1964); A, I. 
Larkin and Yu. N. Ovchinnikov, Sov. Phys.-JETP, 20, 
762 (1965); W. Wu and P. Adams, Phys. Rev. Lett. 73, 
1412 (1994); K. Gloos, R. Modler, H. Schimanski, C. D. 
Bredl, C. Geibel, F. Steglich, A. I. Buzdin, N. Sato and 
T. Komatsubara, Phys. Rev. Lett. 70, 501 (1993). 

[16] F. Braun, Jan von Delft, D. C. Ralph and M. Tinkham, 
Phys. Rev. Lett. 79, 921 (1997). 

[17] Kun Yang and S. L. Sondhi, Phys. Rev. B 57, 8566 
(1998). 

[18] For a mixed phase symmetry that breaks time reversal 
symmetry locally, like d + e l a, where a — Iso-s or d xy 
see last ref. [^| Haranath Ghosh, Europhys. Lett. 43, 707 
(1998); Haranath Ghosh, Phys. Rev. B 59 3357 (1999). 

[19] R. J. Cava, B. B. Batlogg, C. H. Chen, E. A. Rietman, 
S. M. Zahurak and D. Werder, Nature 329 423 (1987). 



7 



